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Abstract. Given an operator L acting on a function space, the J-matrix method consists 
of finding a sequence ?/„ of functions such that the operator L acts tridiagonally on ?/„. Once 
such a tridiagonalization is obtained, a number of characteristics of such an operator L can 
be obtained. In particular, information on eigenvalues and eigenf unctions, bound states, 
spectral decompositions, etc. can be obtained in this way. We review the general set-up, and 
we discuss two examples in detail; the Schrodinger operator with Morse potential and the 
Lame equation. 



1. Introduction 

In many problems one is interested in the eigenfunctions of an operator L acting on some 
function space, or more generally on the spectral decomposition of such an operator in case 
L is moreover self-adjoint. The purpose of the paper is to give an introduction to a method 
that has been used on several occasions and at several places in the literature. This method 
is known as the J-matrix method or as tridiagonalization. A J-matrix, or a Jacobi operator, 
is a tridiagonal operator on some finite dimensional Hilbert space or on the sequence space 
which is usually assumed to be symmetric and having no non-trivial closed reducing 
subspaces. The last conditions are in general not imposed in this paper. A tridiagonalization 
of an operator L acting on some function space is given by a set of functions {?/n}^o ^^^^ 
that L acting on these functions is tridiagonal with respect to these functions, i.e. such that 
(12. ip holds. Note that in the particular case that the upper and lower diagonal term vanishes, 
this just means that the functions ?/„ are eigenfunctions for the operator L. Since there is 
an intimate relation between orthogonal polynomials and three-term recurrence relations, see 
e.g. [E], [26], [29], [36], [38], in such a way that orthogonality properties of the polynomials 
correspond to the spectral properties of the corresponding Jacobi operator, this can then be 
used to find information on eigenfunctions, spectral properties, etc, of the original operator 
L. 

This method is frequently used in physical and chemical models, see e.g. [T], [2], [3], [1], 
[5], [8], [n], [12], [13], [16], [21], [27], [S], [37], [iQ] and references given there. It concerns 
mostly one-dimensional models, and the potentials and Hamiltonians discussed include sextic, 
harmonic oscillator, (Dirac-) Coulomb, (Dirac-) Morse, etc. Usually the papers mentioned 
start out with the operator L to be analyzed, and occasionally with the form of the polynomials 
prescribed, e.g. as in [8] where the Un are monomials times a fixed function. This method is 
also closely related to the Lanczos algorithm in numerical analysis, see e.g. [T^l Ch. 2], and 
to related Krylov subspaces. 
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The purpose of the paper is to discuss the method of tridiagonahzation in a fairly general 
fashion and to consider the special case of a Lame type operator, showing that it can be 
tridiagonalized. This is motivated by the classical theorem of Bochner [TU], recalled in Theo- 
rem \'S.1\ which classifies all orthogonal polynomials that are eigenfunctions to a second order 
differential operator, see also pSl Ch. 20] for generalizations to difference operators, and by 
the classification theorem of Al-Salam and Chihara [7j, recalled in Theorem 13.31 of orthogonal 
polynomials whose derivative can be expressed in a simple way in terms of the orthogonal 
polynomials itself. In general it is difficult to say if an operator can be tridiagonalized, but in 
Section [2] we prove this for a special class of operators including second order differential and 
difference operators with polynomial coefficients of some degree, and we discuss an explicit ex- 
ample in Sections 13.31 and | 3.4[ If there is way to transform, e.g. by conjugation and/or change 
of variables, to such a specific operator, then we can tridiagonalize the resulting operator, as 
is the case for the examples in §21 

It should be noted that a Jacobi operator has simple spectrum, and that conversely a self- 
adjoint operator with simple spectrum can be realized as a Jacobi operator, see [351, Ch. VII], 
assuming that there are no non-trivial (closed) reducing subspaces, see also [9]. This is of 
particular interest in case of the Schrodinger operator, where one can make use of scattering 
theory in order to determine its spectral decomposition. In case both the tridiagonahzation 
procedure works and the spectral decomposition can be made explicit by e.g. an integral 
transformation, the methods can be linked to each other leading to results for the special 
functions and orthogonal polynomials involved and we discuss an example for the Schrodinger 
equation with Morse potential due to Broad and Diestler, see [13], [IE], [TT], [12], [27], as well 
as [30]. 

The contents of the paper are as follows. In Section [2] we discuss a general set-up for 
tridiagonalizable operators. In Section [3] we restrict to second order differential operators, 
where in particular we discuss the Broad-Diestler example and the case of the Lame operator. 

We want to point out that in many cases which are considered there is a link to the bispectral 
problem, see e.g. [2^ for an introduction, and that the tridiagonahzation can be used for both 
the operator in the geometric variable as for the operator in the spectral variable. It is also 
to be pointed out that one can also tridiagonalize (second order) difference operators, which 
are included in the general scheme of Section [2l and that one important example is already to 
be found in Groenevelt [22] for the case of the Wilson functions and the associated difference 
operator. Finally, we want to mention two, closely related, possible extensions that can be 
useful as well. First, one can relate an operator to a doubly infinite Jacobi matrix (i.e. acting 
on £2(Z) instead of on £^(N)), see e.g. [M], [29], and [U Ch. VII]. As indicated by Berezanskii 
[HI Ch. VII] one can also consider this case as 2 x 2-matrix- valued variant of tridiagonahzation, 
and this can then be looked at from a matrix analogue of the tridiagonal situation, see e.g. 
[TH] for an introduction to matrix-valued orthogonal polynomials. This can be useful for such 
operators as the Dirac operator, see [2], [3], [S3] and also [IT] in this context. 

Acknowledgement. Most of this paper was presented by the second author at the meeting 
New Trends in Mathematical Methods for Engineers and he thanks Bruna Germano, Renata 
Martinelli and Clemente Cesarano for organizing this meeting. 
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2. The general set-up 

We consider first a special class of second order operators that can be tridiagonalized, which 
is done in Section 12.11 In Section 12.21 we moreover assume that this operator is symmetric, 
and we consider possible self-adjoint extensions and its spectrum. 

2.1. Motivation and definition. Consider a linear operator L acting on suitable function 
spaces; typically L is a differential operator, or a difference operator. We look for linearly 
independent functions {?/„}5^Lo such that L is tridiagonal with respect to these functions, i.e. 
there exist constants y4„, C„ [n G N) such that 

Lyo = Ayi + BoVq. 

We combine both equations by assuming Cq = 0. It follows that Yl'^=o Pn{z) yn is a formal 
eigenfunction of L for the eigenvalue z if pn satisfies 

(2.2) ZPniz) = Cn+lPn+liz) + BnPn{z) + An-lPn-l{z) 

for n G N with the convention y4_i = 0. In case C„ 7^ for n > 1, we can define po{z) = 1 
and use (12. 2p recursively to find Pn{z) as a polynomials of degree n in z. In case AnCn+i > 0, 
i?„ G M the polynomials p„ are orthogonal with respect to a positive measure on R, and the 
measure and its support then can give information on L in case {?/„}^o gives a basis for the 
function space on which L acts, or for L restricted to the closure of the span {y„}^Q (which 
depends on the function space under consideration). 

We now consider a more specific form of the operator L. Let 5* be a linear operator acting 
on suitable function spaces including the polynomials. We assume that S preserves the space 
of polynomials, and that S lowers the degree by 1, i.e. S = dkX^~^, /c G N, with dk 
for k > 1 and do = 0. Similarly, T is a linear operator acting on suitable function spaces 
including the polynomials. We assume that T preserves the space of polynomials, and that T 
lowers the degree by 2, i.e. T = rf'^x'^"^, G N, with c?'^ 7^ for > 2 and d'^ = d[ = 0. 

Example 2.1. T = 5^ and 5 = ^, the g-derivative S = Dq, or any other g-derivative, see 

e.g. m- 

We now consider the operator L on suitable function spaces 

(2.3) L = MaT + MbS + Mc 

where Mf denotes the operator of multiplication by a function /. We assume that A, B and 
C are fixed polynomials with deg(y4) = a, deg(i?) = h and deg(C) = c. In this case it follows 
that L maps a polynomial of degree n in general to a polynomial of degree max(a + n — 2,h + 
n — l^c + n). So if we look for a tridiagonalization in terms of y„ a polynomial of degree n we 
require a < 3, 6 < 2 and c < 1. 

The case a < 2, 6 < 1 has been studied extensively, in particular the existence of polynomial 
eigenf unctions for Ma T + Mb S for a < 2, 6 < 2, see Bochner's Theorem 13.11 for the classical 
case of T = S"^, S" = and for several other instances of the operators T and see [26l Ch. 
20]. In most of these cases MaT + Mb S have polynomial eigenfunctions which are classes of 
orthogonal polynomials, so that L = Ma T + Mb S + Mq is tridiagonal with respect to these 
polynomials by the three-term recurrence relation in case deg(C) = 1. 
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So the previous discussion motivates why we consider operators as in the following definition. 

Definition 2.2. Let S and T be linear operators preserving the space C[x] of polynomials, 
such that S, respectively T, lowers the degree by 1, respectively 2. We say that the linear 
operator L = Ma T + Mb S + Mc is a TD-operator if A, B and C are polynomials with 
deg(y4) = a < 3, deg{B) = b < 2 and deg(C) = c < 1 with a = 3 or b = 2. Here Mj denotes 
multiplication by the function f . 

Tfieorem 2.3. Let L be TD-operator, then there exists monic polynomials {yn}'^=o, deg(y„) = 
n, such that (12.11) holds for suitable coefficients An, Bn, Cn- 

Proof. First note that there is no loss by assuming the polynomials ?/„ to be monic. 



Recall we assume S x'' = dkx'' ^, k E N, with dk foi k > 1 and do = 0, and T x'^ = 
d'^x''-^, keN, with ^ ^ for A; > 2 and c/q = d[ = 0. Put 



(2.4) Lx'' = {asd'^ + p2dk + 7i) x''^^ + ("24 + + 7o) x'' + {aid'^ + /5o4) ^ + ao4 x^ . 



In particular, the result follows with yn{x) = x" in case ao = 0. 

Now take yo{x) = 1, so that Lyo{x) = C{x). Putting yi{x) = x + co(l) we find that 
Lyo = Aoyi + -BqZ/o if we take Aq = 71, 7iCo(l) + -Bo = 7o- Note that there is choice for 
the constant term Co(l) in yi. Proceeding inductively, we assume that we have determined 
{z/o, ■ ■ ■ ,yk} such that 

(2.5) Lyn = Anyn+1 + Bnyn + Cnyn-l, < U < k - 1, 

Since yk and yk+i are monic polynomials, we see that (12.51) to hold for n = k forces A^. = 
a^d'i. + /?24 + 7i by (12.41) . Putting yk+i{x) = x^~^^ + '^p^^CpX'', we see that we need to 
determine Cp, B^ and Ck from 

(2.6) AkCp = coef[p{Lyk) - Bk coeSp{yk) - Ck coeSp{yk-i), < p < k, 

where coeffp(r) is the coefficient of x^ in a polynomial r. Starting with p = k, we see that 
we need to choose Ck, Bk satisfying — recall yk monic — AkCk = coeSk{Lyk) — Bk, which 
can be easily done for all values of Ak- So we fix Ck and Bk- Next for p = A; — 1 we get 
AkCk-i = coeSk~i{Lyk) — BkCoeSk-iiVk) ~ Ck, for which we choose a solution for Ck-\ and 
Ck- Now we have fixed Bk and Ck, we can solve Cp, < p < k — 2 uniquely (in case Ak 7^ 0) 
from (12. 6p . and we can assign some value to Cp in case Ak = 0. □ 

Remark 2.4. (i) Note that in case e.g. S = and T = S"^ one could stop after the remark 
following (12.41) . since we can use an affine transformation to reduce to the case A(0) = 0. 

(ii) Note that there is freedom in the choice for yn+i in the proof of Theorem 12.31 More 
requirements on the functions should indicate which set to favour. 



A{x) 
B{x) 
C{x) 



03 x'^ + 0:2 + ai a; + ao, 
/52X^ + /?ix + (3q, 

'I1X + 70- 



This implies 
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2.2. Symmetric TD-operators. Now we assume that we have an Hilbert space Ti of func- 
tions containing the polynomials C[x] ^ Ti injectively. We do not assume that L can be 
extended as a bounded operator to Ti, but we assume that L can be viewed as a densely 
defined operator on Ti such that C[x] C D{L), the domain of L. Note that we assume 
L: C[x] — > C[x], and we assume that C[x] dense in Ti. by switching to the closure of C[x] in 
H if necessary. 

Proposition 2.5. Let L be TD-operator. Assume L with domain D{L) is symmetric as 
unbounded operator on Ti, then we can assume {yniUm) = forn ^ m. 

Proof. Since {?/„}^o is a family of polynomials in Ti we can apply the Gram-Schmidt pro- 
cedure to {yo, 1/1, ■ ■ ■ }; and denote the resulting orthogonal set of monic polynomials by 
r„, then we have deg(r„) = deg(?/„) = n and r„ = y„ + Y^kKn^^kUk- By ([2ll]) we find 
Lvn = AnTn+i + Ylik<n^'k''^k- Then {LTn^T^) = for m > n + 1, and for m < n — 1 we 
have 

{Lrn,rm) = {rn, L*rm) = {rn,Lrm) = {rn, Amrm+l + ^ ] ^k^k) = 

k<m 

since m + 1 < n. Note that r.^ G C[a;] C D{L) C D{L*). So L is tridiagonal with respect to 
the orthogonal set {r„}^Q. □ 

Note that Proposition 12.51 easily extends to L skew-symmetric. 

Remark 2.6. Assume that in Proposition 12.51 the orthogonal polynomials yn are eigenfunc- 
tions of a symmetric operator D, Dyn = \iyni such that D preserves the polynomials, 
D: C[x] — i> C[x], and the degree, deg{Dx^) = k, see e.g. Bochner's Theorem 13.11 for the 
classical orthogonal polynomials and, more generally, for all polynomials in the Askey-scheme 
and its g-analogue, see [2H]. So in particular, we assume A„ 7^ 0, n > 1. We assume that D 
acts as a possibly unbounded linear operator on H. Let X be the operator of multiplication 
by the independent variable, so that by orthogonality 

Xyn = an yn+1 + bnyn + C„, yn~l- 

We also assume that X : C[x] — ^ C[x] acts as a possibly unbounded self-adjoint operator on 
7i. Then the anticommutator DX + XD is symmetric, and by 

{DX + XD) yn = a„(A„+i + A„) y^+i + 2A„6„ y„ + c„(A„ + A„_i) 

it follows that DX + XD is a symmetric TD-operator. 

Conversely, if L is as in Proposition 12. 5[ then we can define D as a linear operator on C[x] 

by 

n-l 

(2.7) Dx" = ^(-l)*^X*^Lx"-^-^ n > 1, D1 = 0, 

k=0 

by iterating D x^ = DX x"'~^ = (L — XD)x"'~^ and using the initial condition D 1 = 0. Note 
that this completely determines D on the polynomials C[x]. From fl2.7p we can show that 
DX + XD = L on C[a:]. Since we assume L and X symmetric, we get D*X + XD* = L on 
C[x] assuming D* preserves the polynomials. If one also assumes that degD*x'' < k, we see 
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that D* must have the same form as D. So D is symmetric if we can show that D*l = 0. By 
the assumptions we have D*l = c for some constant c, and 

_ {D*l,l) _ {l,Dl) 

Since D is symmetric, preserving polynomials and the degree, we find Dyn = XnUn for real 
\n with Ao = 0. 

In case L is antisymmetric, this has been completely worked out by Koornwinder [311 §2], 
and then one has interesting links to the so-called string equation. 

In the situation of Proposition 12.51 we can next ort honor malize the orthogonal polynomials 
{?/„}.^o then we get 

(2.8) Lyn = An Vn+i + Bn y„ + A^-i yn-i 

with An, Bn G M and with the convention A^i = 0. Note that in the skew-symmetric case we 
obtain the same result but with An, Bn G iM and with the convention A^i = 0. 

The situation in (12. 8p is governed by the occurrences of An = 0. In case An-^^ = and 
An2 = with ni < n2, and, in view of the convention, ni = — 1 is allowed, we see that L 
preserves the finite-dimensional subspace spanned by ?/„ for ni < n <n2, which has dimension 
ri2 — ni . In particular, if n2 = ni + 1 we see that yn2 is an eigenfunction of L for the eigenvalue 
Bn2- We have to distinguish between the cases of finite or infinite zeros oi n ^ An- 

Theorem 2.7. Let {L,D{L)), with D{L) = C[x] be a symmetric densely defined TD- 

operator with the tridiagonalization (12. Sp . Assume — 1 = no < < n2 < • • ■ is such that 

An, = 0, 

(i) In case N 3 n An has an infinite number of zeros, the finite- dimensional subspaces 
Tii = span{?/„ | Ui^i < n < Ui}, i > 1, dimTij = Ui ~ Ui^i, are invariant for L. Moreover, 
Ti = Til and L\-ii^. has simple spectrum consisting o/dimTij different eigenvalues. The 
operator {L,D{L)) is essentially self-adjoint. 

(ii) In case N 3 n \—>- An has k zeros, = —1 < ni < ■ ■ ■ < n^, the k finite dimensional 
subspaces Hi = span{|/„ | ?T,j„i < n < n^}, 1 < i < k, dimTij = Ui — Ui^i, are invariant for L. 
Llfi- has simple spectrum consisting o/dimTij different eigenvalues. Consider the sequence of 
polynomials determined by pq{z) = 1, and 

ZPn{z) = An+nk+lPn+l{z) + Bn+nk+lPn{z) + A„+„^p„_i(z) 

then {L,D{L)) is essentially self-adjoint if and only if the orthogonal polynomials {pn}'^=o 
correspond to a determinate moment problem. 

Proof. In case = and An^ = with ni < n2 we see that L preserves the finite-dimensional 
subspace /C, dim/C = n2 — ni, spanned by ?/„ for ni < n < n2- By (12. Sp it follows that 
L : K, ^ K, is given by a Jacobi matrix, i.e. a symmetric tridiagonal matrix. It is well-known, 
see e.g. that such a matrix has dim/C different eigenvalues, and that each of them 

has multiplicity one. In case (i) we have that the closure of L is given by its maximal extension, 
which is self-adjoint. 

In case (ii) the previous considerations remain valid for the finite dimensional invariant 
subspaces, and we are left with the study of the action of L on the closure JC of the hnear 
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span {yn+nk}'^=o- Let £^(N) be the Hilbert space of square summable sequences with standard 
orthonormal basis {e„}^Q. Then U : fC ^ £^(N), Vn+rn, i— >■ e„ is a unitary map such that 

So the action of L restricted to /C is intertwined with the action of a Jacobi operator on £^(N), 
and it is well-known, see e.g. [22], [35j, ^38j, that this Jacobi operator is essentially self-adjoint 
if and only if the corresponding moment problem is determinate. □ 

The spectrum of a TD-operator on finite-dimensional invariant subspaces can be determined 
explicitly, and in case we can also find the eigenfunctions in another (direct) way this leads 
to non-trivial sums, see e.g. §3.31 for an example. Let us now assume that the TD-operator 
L with domain D{L) acting on Ti. with D{L) = C[x] ^ Ti dense in 7^ = L^(z/) is essentially 
self-adjoint, and we assume that A„ has no zeros, except the convention A_i = 0. So we are 
in the second case of Theorem 12. 7[ In this case the spectrum is simple [33, Ch. VI] , so the 
spectral theorem states that there exists a unitary map T : Ti = [v) ^ K, = -L^(yu), to some 
weighted L^-space with a positive Borel measure on R such that T L T* = X, where X is 
the (possibly) unbounded operator on L'^ifJ^) of multiplication by the independent variable, say 
A, see |35l Ch. VI]. We assume that there exist suitable functions 0a, generally not assumed 
to be in the Hilbert space H, such that (T/)(A) = {f,(f)x) for suitable f G H and where 
L(j)x = A0A- This is a typical situation in the spectral decomposition of various second order 
differential or difference operators. 

In this case L has simple spectrum, and since T?/„ satisfies the same recurrence relation we 
find 

(2.9) (Ty„)(A) = J ct>x{x)yn{x)dv{x) = p„(A) (T1)(A), 

or, the integral transform with kernel the (formal) eigenfunctions of L maps the orthogonal 
polynomials ?/„ to the orthogonal polynomials p„, up to a common multiple. See e.g. [22], [SD] 
for examples. 

3. Second order differential operators 

We now restrict ourselves to the case of second order differential operators as an example. 
Needless to say that appropriate g-analogues or difference analogues can be considered as well 
within this general framework. First we discuss some generalities, and then we discuss two 
examples; the Schrodinger equation with the Morse potential in Section 13.31 and the Lame 
equation in Section [3^ 

3.1. Theorems by Bochner and Al-Salam— Chihara. We now assume that 
(3.1) L = Ma^ + Mb^ + Mc, 

so we take S = and T = S"^ = This then fits into the scheme of Theorem 12.31 Recall 
our basic assumption that deg{A) = a < 3, deg(-B) = b < 2, and deg(C) = c < 1, and that 
we assume that a = 3 or 6 = 2. Indeed, in case a < 2 and 6 < 1 we are essentially back 
to Bochner's Theorem 13.11 and the fact that all polynomials in Bochner's Theorem satisfy a 
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three-term recurrence. For completeness, we recall Bochner's Theorem here, see Bochner |lUj . 
or e.g. §20.1]. 

Theorem 3.1 (Bochner (1929)). Up to affine scaling the only sets {?/„}^q of polynomials that 
are eigenfunctions to a second order differential operator A{x) y'nix) + B{x) y'n{x) + A„ yn{x) = 
are 

(1) Jacobi polynomials: deg(y4) = 2 with 2 zeroes, deg{B) = or 1; 

(2) Laguerre polynomials: deg(y4) = 1, deg(i?) = 1; 

(3) Hermite polynomials: deg(A) = 0, deg(i?) = 1; 

(4) Bessel polynomials: deg(A) = 2 with double zero, deg{B) = or 1 and A and B have 
no common zero; 

(5) Monomials: deg(y4) = 2 with double zero, deg(-B) = 1 and A and B have a common 
zero. 

For the notation of orthogonal polynomials we follow the notation as in [28], and for the 
Bessel polynomials we follow [26j, so P!f''^\x), Li^\x) and Hn{x) denote Jacobi, Laguerre 
and Hermite polynomials, whereas yn{x;a,b) are Bessel polynomials. In Bochner's Theorem 
13.11 the first three sets are orthogonal polynomials on the real line and these are classical 
orthogonal polynomials. The Bessel polynomials are not orthogonal on the real line with 
respect to a positive measure, see [26j §4.10], and the same is true for the monomials. 

Remark 3.2. Bochner's theorem has several analogues, e.g. by replacing the differential 
operator ^ by one of the g-difference operators, see [261 Ch. 20] for more information. 

Next one can ask for a relation between the derivative of an orthogonal polynomial, and its 
relation to orthogonal polynomials of possibly different degree within the same family. The 
following classification theorem has been obtained by Al-Salam and Chihara [7], see also the 
survey by Al-Salam [6]. 

Theorem 3.3 (Al-Salam and Chihara (1972)). If {pn}^=o is a family of orthogonal polyno- 
mials on the real line with differential-recursion relation 

Gix) ^(x) = AnPn+l{x) +BnPn{x) +CnPn~l{x) 

ax 

for constants An, Bn, Cn and G (necessarily) a polynomial of degree < 2, then the Pn's are 
(up to affine scaling) Jacobi, Laguerre or Hermite polynomials: 

(1 - x') j-/'n^'\x) = Al^'^^ Pi:f (x) + P(-^)(X) + C("'^) Pi-f)(x); 

= nLl-\x) -in + a)Li"l,ix); 

-^Hnix) = 2nHn-i{x). 
ax 

Remark 3.4. The classification of Al-Salam and Chihare concerns orthogonal polynomials, 
so that the Bessel polynomials and the monomials do not occur in the list. However, the 
Bessel polynomials and the monomials satisfy a differential-recursion relation of the form 

x^^yn{x]a,b) = A'^'' yn+iix; a,b) + B^''' ynix; a,b) + C^'^ y^-iix; a,b), -^x" = nx"-\ 
ax ax 
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So Bochner's Theorem 13.11 and Al-Salam's and Chihara's Theorem 13.31 deal with the same sets 
of polynomials. 

As noted in Remark I2.4( i). by an affine transformation we can assume ^(0) = in fl3.ip . 
and then L is tridiagonalized by the monomials, cf. Theorem 12.31 and its proof. However, 
there is choice in the polynomials leading to tridiagonalization. Using Bochner's Theorem 13. II 
and Theorem 13.31 by Al-Salam and Chihara and Remark 13.41 one can proceed as follows to 
tridiagonalize the operator L: first use Bochner's Theorem 13.11 to get rid of the second order 
derivative; second use Theorem 13.31 to get rid of the the first order derivative. Note there is 
a lot of choice: first by using an affine transformation; second by choosing the parameters in 
case of the Jacobi, Laguerre, or Bessel polynomials; and thirdly in the possible decomposition 
of the polynomial A as a product of two lower order polynomials. We give an example of this 
procedure when discussing the Lame equation. 

3.2. Symmetric second order differential equations. We now consider the case of L = 
Ma-^ + Mb + Mc being symmetric on a Hilbert space H = L'^(^{a,b),w{x)dx), where 
— oo < a < 6 < oo and w{x) > on (a, 6). Recall that we assume C[x] "-^ Ti as a dense 
subspace. 

Lemma 3.5. Assume A, B and C are real-valued polynomials on R. Moreover, assume 
w G C^{a,b), and (Aw)' = Bw, then L, with domain D{L) consisting of C^{a,b), is a 
symmetric operator onTC = L'^[{a, b), w{x)dx^ . 

Proof. For f,g& (a, b) we have 



{Lf){x)g{x)w{x)dx= I [C{x)f{x) + B{x)f'{x) + A{x)f"{x))g{x)w{x)dx 



(3.2) 



C{x)f{x)g{x)w{x)dx+ / f'{x)B{x)g{x)w{x)dx 

J a 

pb pb 

— / A{x)f'{x)g'{x)w{x)dx— / f'{x)g{x)(^Aw) (x) dx 



C{x)f{x)g{x)w{x)dx— / A{x)f'{x)g'{x)w{x)dx 

J a 

since [Aw)' = Bw. The right hand side of (13. 2p yields the symmetry. □ 

Assuming the conditions of Lemma [3.51 on the weight function w, we can write, for f,g& 
C[x], 



{Lf,g)= / C{x)f{x)g{x)w{x)dx— / A{x)f'{x)g'{x)w{x)dx 

J a J a 

+ A{b)wib) f{b)g{b) - A{a)w{a) f{a)g{a), 

so that Lemma ES] has the following analogue in case the domain D{L) = C[x] is considered, 
in case 
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Lemma 3.6. Assuming the conditions of Lemma \3.5\ on A, B, C andw. Moreover, assume 
Aw has a zero in a and b, which has to be interpreted for a = oo, respectively b = —oo, as 
\imx^ooU!{x)p{x) = 0, respectively \\mx-*~oow{x)p{x) = 0, for all polynomials p . ThenL, with 
domain D{L) = C[x], is a symmetric operator onTi. = L'^[{a,b),w{x)dx^ . 

The first order differential equation for tlie weight function w in Lemma 13.51 is rewritten as 

(3.3) {\nw) =- = — 

i.e. the logarithmic derivative of is a rational function for which the degree of the nu- 
merator polynomial is at most 2 and the degree of the denominator polynomial is at most 
3. Depending on the structure of the rational function the differential equation (13. 3p can be 
solved straightforwardly using a partial fraction decomposition. The solution of (13. 3p will very 
much depend on the relation between the polynomials A and B. 

E.g. in the special case A' = B we see that we can take w{x) = 1, and Lemma applies, 
and from Lemma [3.61 we see that L with D{L) = C[x] is symmetric on L'^[{a, b), dx) if a and 
b are different zeroes of the polynomial A, a < b. 

3.3. Schrodinger equation with Morse potential. The Schrodinger equation with Morse 
potential is studied by Broad [13j and Diestler [16j in the study of a larger system of cou- 
pled equations used in modelling atomic dissocation. The Schrodinger equation with Morse 
potential is used to model a two-atom molecule in this larger system. 
The Schrodinger equation with Morse potential is 

(3.4) -^ + ^' qix) = b\e-'^-2e-n, 

which is an unbounded operator on L^(R). Here 6 > is a constant. It is a self-adjoint operator 
with respect to its form domain, see [HH Ch. 5] and lim2;^oo = 0, and lim^^^-oo q{x) = +oo. 
Note min(g) = —6^, so that discrete spectrum is contained in [—6^, 0] and it consists of isolated 
points. We look for solutions to —f"{x) + <i{x)f{x) = 7^/(x). Put z = 2be~^ so that x G M 
corresponds to z G (0, oo), and let f{z) = -^gi^z), 

(3.5) + (Zif!±^i±l!±i),(,) . 0, 

Z'^ 

which is precisely the Whittaker equation with k = b, fi = ±Z7, and the Whittaker integral 
transform gives the spectral decomposition for this Schrodinger equation, see [IHl § IV]. In 
particular, depending on the value of b the Schrodinger equation has finite discrete spectrum, 
i.e. bound states, see the Plancherel formula fl9[ § IV], and in this case the Whittaker function 
terminates and can be written as a Laguerre polynomial of type Lm^~'^™'~^\x) , for those m G N 
such that 2b — 2m > 0. 

The Schrodinger operator is transformed into a TD-operator, and a particularly nice basis in 
which the operator is tridiagonal is obtained by Broad [13] and Diestler (TB]. Put N = #{n G 
N I n < b— |}, i.e. = [b+ |J, so that 2b — 2N > —1, and we assume for simplicity b ^ | + 

Let T: L'^{R) L^{{0,oo);z^^~^^e~'dz) be the map (T/)(z) = z^-^-h^^' f{\n{2b/z)), then 



T is unitary, and T(-^ + q)T* = L with L = Maj^ + Mg^ + Mc with A{z 



z 
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B{z) = {2N-2b-2 + z)z, C{z) = -{N-b-lf + z{l-N). Using the second-order differential 
equation, see e.g. [26], (4.6.15)], [28], (1.11.5)], [Ml (5.1.2)], for the Laguerre polynomials, cf. 
Bochner's Theorem 13.11 the three-term recurrence relation for the Laguerre polynomials, see 
e.g. [261 (4.6.26)], [28, (1.11.3)], [36, (5.1.10)], and the differential-recursion formula as in 
Theorem 13.31 for the Laguerre polynomials we find that this operator is tridiagonalized by the 
Laguerre polynomials Ln'^~'^^\ When we translate this back to the Schrodinger operator we 
have started with we obtain 



(3.6) yn{x) = (26)(^-^+i) J— -{h-N+\).^-be-^ Lf -2^)(26e-^) 

^ ' yn\ ) \ ) Y r(26-2Ar + n + 1) " ^ ' 

as an orthonormal basis for L^(]R) such that 

+ ijVn = -{1-N + n)^{n + l){2h-2N + n + 1) y^+i 

^^■^^ + (^-(iV - h - ]^f + {I - N + n){2n + 2h - 2N + I) - yn 

-in- N)^n{2b-2N + n) y^-i. 

Note that (13.71) is written in a symmetric tridiagonal form. 

The space spanned by and the space H~ spanned by {yn}n=d invariant 

with respect to + g which follows from Note that ^^(R) = H+^H-, dim(7^~) = N. 

In order to determine the spectral properties of the Schrodinger operator in this way we 
follow the approach of Theorem 12.71 We first consider its restriction on the finite-dimensional 
invariant subspace 7i~. We look for eigenfunctions "^^Iq Pn{z) yn for eigenvalue z^ so we need 
to solve 



zPn{z) ={N -l-n)^{n + l){2h -2N + n + l) Pn+i{z) 

+ (-(iV - 6 - i)2 + (1 - AT + n)(2n + 26 - 2Ar + 1) - Pn{z) 

+ {N - n)^n{2h -2N + n) P„_i(z), < n < iV - 1. 

which corresponds to some orthogonal polynomials on a finite discrete set. These polynomials 
are expressible in terms of the dual Hahn polynomials, see §6.2], [^ §1-6], and we find 
that z is of the form —{h — m — |)^, m a nonnegative integer less than 6 — |, and 



P„(-(6 - m - 1)2) = ^ ^J^^^" Rn{X{N - 1 - m); 26 - 2iV, 0, AT - 1), 

using the notation [26l §6.2], [28], §1.6]. Since we have now two expressions for the eigenfunc- 
tions of the Schrodinger operator for a specific simple eigenvalue, we obtain, after simplifica- 
tions, 

N-l 

Rn{X{N-l-m);2b-2N,0,N-l)L^^'-^''\z) = C z""-^-^ L^^'-^'^-^\z), 

(3.8) 

C = (-1)^+-+! ({N + m- 2b)N-i-m ~ ^ 
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where the constant C can be determined by e.g. considering leading coefficients on both 
sides. Using the orthogonahty relations [2H1 (1-6.2)] of the dual Hahn polynomials, 03.81) can 
be inverted. 

On the invariant subspace we look for formal eigenvectors of the form Yl^=o Pn{z) yN+n{x) 
for the eigenvalue z. This leads to the recurrence relation 



Pniz) = -(1 + n)^/iN + n + 1)(26 - N + n + 1) Pn+i{z) 

+ {-{N - 6 - i)2 + (1 + n)(2n + 26 + 1) - n - A^) P„(z) 



n 



v/(iV + n)(26-iV + n)P„_i(z). 



This corresponds with the three-term recurrence relation for the continuous dual Hahn poly- 
nomials, see [2Sj §1.3], with (a, 6, c) replaced by {h + \, N — h + — N + ^) , and note that 
coefficients a, h and c are strictly positive. We find, with 2; = 7^ > 



Pn{z) 



n\^{N + l)^{2h-N + l) 
and 



00 

2\ D 1^.2 



Pnil )Pmil ) ^(7) d-f = Sn,m, 

r(6 + i + zj)T{N - 6 + i + «7)r(& - iV + i + 27) 



27rA^!r(26- 1) 



r(227) 



Note that the series Yl'^=o Pnil^) yN+n diverges in 7i+ (as a closed subspace of L^(R)). Using 
the results on spectral decomposition of Jacobi operators, we obtain the spectral decomposi- 
tion of the Schrodinger operator restricted to Ti.'^ as 

T-.n^^ L2((0,oo);^(7)rf7), (T|/,v+n)(7) = Pnil'), 

{{-^,+q)f,9) = I 7'(T/)(7)(T^7)(7)^(7)^7 

for fiQE Ti.^ C L^(]R) such that / is in the domain of the Schrodinger operator. 

In this way we have obtained the spectral decomposition of the Schrodinger operator on 
the invariant subspaces H~ and H^, where the space H~ is spanned by the bound states, 
i.e. by the eigenf unctions for the negative eigenvalues, and 7i+ is the reducing subspace on 
which the Schrodinger operator has spectrum [0, 00). The link between the two approaches 
for the discrete spectrum is given by (13.81) . For the continuous spectrum it leads to the fact 
that the Whittaker integral transform maps Laguerre polynomials to continuous dual Hahn 
polynomials, and we can interpret (13. 8p also in this way. For explicit formulas we refer to [SOI 
(5.14)]. Koornwinder [30] generalizes this to the case of the Jacobi function transform mapping 
Jacobi polynomials to Wilson polynomials, which in turn has been generalized by Groenevelt 
|22] to the Wilson function transform mapping Wilson poljTiomials to Wilson polynomials. 
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3.4. Lame equation. The classical Lame equation is — (m(m+ l)p{u) + E^F{u) = 0. 

Here p is the WeierstraB p-function, which is a doubly-periodic function with periods 2ui, 
2uj2 (and ^ ^ M). We do not yet assume a condition on m e M, but note the symmetry 
m —m — 1. This equation is very classical, and it is studied in |39« §23] in detail. Put 

X = p{u), and F{u) = f{p{u)) then 

(ff df 1 

A{x)-^{x) + B{x)-^{x) - -{m{m + l)x + E)f{x) = 0, 

(3.10) A(x) = (x - ei)(x - e2)(x - es), 

B{x) = \{(^~ e2)(a; - 63) + (x - ei)(x - 63) + (x - ei)(x - 62)^. 

Note that the e^'s are all different, where we follow the notation as in Whittaker and Watson 
[39i §20, §20.32] for p, e^, etc. In the form fIXTUD it is a TD-operator. In [39i §23.41] a related 
procedure is discussed which leads to solutions of fl3.10p for specific values of E by inserting 
descending power series in x — 62 . Another classical line of study is to allow for a degree p+1 
polynomial in front of the second order derivative and a degree p polynomial in front of the 
first order derivative in fl3.10p and next to look for a polynomial, known as the Van Vleck 
polynomial, of degree p — 1 in front /(x) in f l3.10p such that f l3.10p has a polynomial solution 
S'(x), known as the Heine-Stieltjes polynomial, see [361 §6-8] ^ind references for this line of 
considerations. 

As noted (13.101) is a TD-operator, which we now tridiagonalize. In light of Bochner's 
Theorem 13. H the Al-Salam and Chihara Theorem 13.31 and the procedure sketched in §3.11 we 
first use an affine transformation x = ay + b, a = i(ei — 62), b = ^{ci + 62), so that y = I 
corresponds x = ei, y = —1 corresponds x = 62. Note that we can use any other permutation 
of the points ei, 62 and 63. This yields 

(y -l)iy+ l)iy - + l{(y + ^)(y - a) + (y - l)iy - a) + (y - l){y + 1)) ^(y) 

- ^{m{m + l){y+^)+ E)g{y) = 
with a = _a+£2ii2e3 ^ _3e3_ / f \ ^ f{—). Let us denote the second order 

61—62 61—62 ' ' J \a I \ a ' 

differential operator for = by L. In view of Bochner's Theorem 13.11 and the factor y^ — \ 
in front of the second order derivative, we try to tridiagonalize the operator using the Jacobi 
polynomials Pi"'^'* and its second order differential equation, see p6], (4.2.6)], [281 (1-8.5)]. In 
this way we get rid of the second order derivative, and collecting the remaining terms in front 
of the first order derivative gives 

{y - a) ((/? - a) - y(« + /5 + 2)) + i ((y + l)(y - «) + (y - _ «) + (y - l)(y + 1)) , 

so that we can use the Al-Salam and Chihara Theorem 13. 31 in case this is a multiple of (y^ — 1). 
So this expression has to be zero for y = ±1, and we find a = /3 = — ^, i.e. we have to take 
the Chebychev polynomials T„ in order to tridiagonalize the Lame equation. So using the Al- 
Salam and Chihara Theorem 13.31 and the three-term recurrence for the Chebychev polynomials 
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Tn, see e.g. fl^ §4.5], [2HI (1-8.34)], we obtain 

LTn= - {2n — m){2n + m + l)Tn+i + (—an^ — -m{m + 1)-')T„ 

(3.11) ^ ^ V 4 a/ 

+ -(2n + m){2n - m - l)^^^! 
o 

for n > 1 and for = 

(3.12) LTo = -i(m(m + 1)(|/ + -))To = -i(m(m + l))Ti - -im{m + l)-)ro. 

4 a 4 ^ 4 ^ a 

Note that we cannot consider 03.121) as the special case n = of fl3.1ip . Note also that 03.111) 
and 03.121) exhibit the symmetry m —m — 1. It is to be noted that the recurrence 03.111) 
can be solved using the continuous dual g-Hahn polynomials, see [28l §1.3], precisely for the 
excluded(!) values a = ±1. Now for the Lame equation we need to solve Lip = E ip. 

Remark 3.7. It should be noted that this relation between the Lame operator and the 
Chebychev polynomials is conceptually different from a link discussed in Finkel et al. |20j . 
which is related to the results of Ince [25]. Finkel et al. [20] use the Jacobian version of 
the Lame operator, whereas we use the algebraic form, see [39l §23.4]. Their approach is 
motivated from the theory of quasi-exactly solvable Hamiltonians, see [21] for an overview. 

From 03.111) it is clear that the Jacobi matrix for the Lame operator splits in case m G Z. We 
only discuss the case m G N is even, since we obtain a finite dimensional invariant subspace. 
This is done in Section [3.4. 1[ In Section [3.4.2l we consider a special case in which no coefficients 
in 03. lip . 03.120 vanish and such that we can write L in a symmetric form. 

3.4.1. Case m = 2k & N is even. Let us first consider the case oi m = 2k is even, then the 
Lame operator L leaves the k + 1-dimensional space spanned by T„, n = 0, . . . , fc, invariant. 
We can rewrite 03. lip in this 

LT„ =(^(2n + l)(2r2) - ^{2k)i2k + 1))t„+i + [-an^ - hk{2k + 1)^)t„ 

(3.13) + Q(2n)(2n-l)-^2A;(2A; + l))T„_i, n > 1, 

1 I h 

= - -(m(m + l))Ti(?/) - -(m(m + l)-)To(y), n = 0. 

We now look for eigenfunctions L ^^=o Pn{E)Tn = E Ylt=o Pn{E)Tn, so we need 

EPo{E) - Uk{2k + l)y,{E) - ^(m(m + l)^)n(^), 

EP^{E) =(^(2n + 2){2n + 1) - Uk{2k + l))p„+i(^) + [-an" - -2k{2k + l)^)^n(^) 

+ i^{2n - l){2n - 2) - ^(2A;)(2A; + l))p„_i(E), l<n<k, 
EPk(E) =(^-ae - -2k{2k + l)^)^fc(^) + (^(2A; - 1)(2A; - 2) - ^(2A;)(2A; + l))p,_i(E). 



THE J-MATRIX METHOD 



15 



The possible values for the eigenvalue E are determined as follows; generate the polynomials 
Pn by the first two equations starting with Po{E) = 1. Stop at Pk+i{E), and then the zeroes 
of Pk+i{E) are the only possible eigenvalues of the Lame operator restricted to this finite- 
dimensional space. This corresponds nicely to §23.41]. Note moreover that for a G M 
Favard's theorem is valid, implying that the polynomials Pn, n = 0, . . . , k, are orthogonal with 
respect to a measure on the real line, so that there are k + 1 different real eigenvalues E. 

3.4.2. Orthonormal version. Assuming that there are only non-zero coefficients in the three- 
term relation (13.111) . (13.121) we can ask under what conditions there exists an orthonormal 
version. Assume m G {2k + 1, 2A; + 2) for some G N, and put T„ = a^Pn with 

/ (|(1 -m))„(l + \m)n 
"""V (-lm)„(i(m + l))„ ■ 

Note that the condition on m implies that the argument of the square root is indeed positive. 
Then (13.111) . (I3.12p is rewritten as 

Lpn=an Pn+l + KVn + On-l Pn-1 



1 L 1 ,,1 1,, 1 1 L 

(3.14) an=-y{n + -m+ l){n - + -){n - -m){n - + -), 

2 1 , ,b 

— an mim + 1) — , 

4 ^ ^ a 

Lpo = 2aoPi + bopo, 

which can be viewed as a symmetric operator assuming a G M except for the the last line in 
(I3.14p . At this point it is not clear if the Jacobi form of L as displayed by the first equality in 
(I3.14P gives rise to an essentially self-adjoint operator or not, since the coefficients a„ = 0{n'^), 
bn = C(n^) blow up. 

Note that in this case (13.31) gives 

B-A' 1/1 1 1 



(Inw)' , I . . 

A 2\y -I y + l y - a 

so that we take if to be a multiple of 



v/(y2 - 1) 

which we might take on [—1, 1] and assuming that y — a is positive on [—1, 1], or a G M\ [—1, 1]. 
A further study of the orthonormal case seems to be required. 
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